Effects of ion and electron correlations on neutrino scattering in the infall phase of a 

supernova. 
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Many authors have used one-component plasma simulations in discussing the role of ion-ion corre- 
lations in reducing neutrino opacities during the collapse phase of a supernova. In a multicomponent 
plasma in which constituent ions have even a small range of N/Z ratios neutrino opacities are much 
larger, in some regions of parameters, than for the case of a one component plasma. 
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There exists an extensive literature concerning the role 
of ion-ion correlations in reducing the neutrino opacity 
in the region that collapses to form the supernova core 
0-0- These correlations are important, by virtue of 
the fact that typical neutrino wavelengths are large com- 
pared to the Debye length. Beginning with the stan- 
dard model result for the coherent neutrino scattering 
from a single ion of charge Z and neutron number N, 1 



part of the neutral current couplings to neutrinos, as, 



where 



d<jQ /dCl 



C = -2Zsin 2 6L 



AO/2, 



(1) 



(2) 



it is standard to express the ionic correlation effects 
through a structure function S(q). Then we write a dif- 
ferential scattering rate for the neutrino, dj/dCl, as, 



d r y{q)/dVl — (dao/dCl)m S(q) , 



(3) 



where nj is the ion number density and q is the momen- 
tum transfer to the ions. 

The calculations of S(q) reported in Q-0 are all based 
on a one component ionic plasma, the electrons being suf- 
ficiently degenerate in the regions of interest as to form 
a virtually uniform background. In ref.Q the authors 
do attempt to simulate the multicomponent case by us- 
ing effective averaged parameters in the one-component 
model. But the present note provides evidence that in 
the regions in which S(q) << 1 and in which the nuclei 
have significant diversity in their N/Z ratios, these one- 
component plasma simulations give neutrino scattering 
rates that are much smaller than those that would be 
obtained in a genuine multicomponent calculation. 

We illustrate this for the case of a two component 
plasma with average number densities n%, n 2 , and nuclear 
charges Z\, Z 2 . Taking into account the weak currents of 
the ions alone, we write the operative vector (coherent) 



Hi = G 



w 



J d 3 xipl(x)ip lJ (x){\ 1 Z 1 n 1 (x) + X 2 Z 2 n 2 (x)\ , 



(4) 



where 



Ai, 2 = — [-2Zi, 2 sin 2 Q w + (Zi j2 - N lt2 )/2] , (5) 

^1,2 

and ri\ (x) , n 2 {x) are the respective density operators for 
the two varieties of ions. We obtain, 



dj(q)/dn 



G 2 w E 2 (l+cost 
4tt 2 



l 3 x e iq x 



Z 2 A 2 (ni(x)ni(0)) + 2X 1 X 2 Z 1 Z 2 (n 1 (x)n 2 (0)) 



-X 2 2 Z 2 {n 2 (x)n 2 (0)) 



(0) 



1 We take h, c, kg = 1 in what follows. 



The fact that Ai is different from A2 for different nu- 
clear species, even if only slightly different, is the key 
to what follows. The physical point that neutrinos scat- 
ter from the fluctuations of the quantity on the RHS 
of (@J, [XiZiriiix) + X 2 Z 2 n 2 (x)}. If Ai = A 2 then the 
fluctuations of this source strength are proportional to 
the fluctuations in electric charge density, and the charge 
density does not like to fluctuate in the q — 0, or long 
range, limit. ( For the moment we assume that the high 
electron degeneracy prevents any fluctuation of the elec- 
tron density.) In the case of the two component plasma 
with Ai 7^ A 2 , the source strength for neutrino scatter- 
ing, which is not proportional to the charge density, can 
fluctuate while leaving the ionic charge density strictly 
zero. 

To address this analytically, we begin in the q — > 
limit, where a simple argument based on statistical me- 
chanics suffices. Then we shall give the solution of the 
two component problem for all q at the Debye-Hiickel 
(D-H) level of approximation. For orientation, we begin 
with the statement from ref.Q that the "correct small q 
behavior" of the structure function S is given for the one 
component plasma by, 



2 



S(q) 



sr 



(aiq) 2 kBTxdn/T 



= l- + 0(q*), (7) 



where T is the conventional plasma coupling constant and 
ai is the mean interionic spacing. In the final equality 
we have substituted the expression for T in terms the 
Debye wave number n 2 = e 2 Z 2 ni/T. The limit J7J illus- 
trates the reluctance of the single component plasma to 
fluctuate. 

For the two-component plasma we define the partial 
squared Debye wave numbers for the respective species as 
Ki 2 = e 2 Zf 2 n\ i/T. In this case (J7J) should be replaced 
by, 



lim S(q 2 

q^O 



T(Xi 



A 2 ) 2 



2„2 



l rv 2 



1 (k 2 + K 2 )mC 2 



(8) 



where C can now be taken as any average coupling con- 
stant factor for the two species. Note that the factor 
ii] C 2 in the denominator cancels when we calculate the 
rate from 

Eq. ||SJ| can be derived from the basic principles of sta- 
tistical mechanics, following the steps of section 2 , "the- 
ory of multicomponent fluctuations" , of ref. ,8§ and en- 
forcing in addition a constraint of complete local neutral- 
ity of the plasma. We only sketch those considerations 
here, beginning from the construction of the limit, in a 
purely classical treatment, of the mean of quadratic forms 
in the Fourier components of fluctuations 5rii(q), 



lim {5ni(q)6rij(- 



•<?)} 



T^(VoL). 



(9) 



For the free energy density functional needed to eval- 
uate the derivatives with respect to the chemical poten- 
tials, fii, we take just the kinetic term plus the simplest 
term that ensures complete local neutrality when a pa- 
rameter, b — > oo. This interaction is an energy per vol- 
ume of V = b(Z\n\ + Z 2 n 2 — n e ) 2 /2. The densities of 
species #'s 1 and 2, are now given by, 



ni — 2 exp 



T 1 [fj,i - bZi(Zini + Z 2 n 2 - n e ) 
d 3 k 



■exp(-fc 2 /2MT), (10) 



and 



ri2 = 2 exp 



T 1 [fj, 2 - bZ 2 {Z 1 ni + Z 2 n 2 - n e )] 
d 3 k 



■exp(-fc J /2MT). (11) 



After taking the logarithms of (|10fl and l|ll|) we dif- 
ferentiate to obtain the matrix dni/drij, which is then 
inverted to provide the RHS of 10, which when substi- 
tuted in © yields the result ©in the b — > oo limit. 



To address the q 2 dependence at the D-H level we begin 
by defining a set of density correlators, 

K id {q) = T- 1 / cfr^-MxJn^O)) , (12) 

as the polarization functions, where T is the tempera- 
ture. In the absence of the Coulomb interactions of the 
particles in the plasma we would have simply Kij = 
T 6i t j(rii}, that is, a diagonal matrix with number den- 
sity values on the diagonal. 

It is almost always useful to define "proper" polar- 
ization parts; graphically in a perturbation development 
for the correlator these are sums of graphs for Kij that 
are individually not divisible into unconnected pieces by 
cutting one Coulomb line. The reconstitution of the 
complete correlator from these proper polarization parts, 
which we designate as 11, j, is fairly straightforward, and 
there are various languages for carrying it out. A com- 
plete and elegant derivation and statement of the results, 
in the case of a classical plasma, is given by Brown and 
Yaffe 0. We quote eq.(2.110) of this paper, with minor 
changes of notation, for the case of any number of ionic 
species, 2 



e a ebK ab (q) = e a e 6 II ab (q) - 
[ e«e c n ac (g)] [ e b e c Ii b M] G(q) , (13) 



where, 



G- 1 (q)=q 2 + J2^e b U ab (q), (14) 

a , b 

and e = eZ a is the charge of nuclear species a. To ob- 
tain the Debye-Hiickel level result we insert the lowest 
order answer for the proper polarization parts and eval- 
uate them at q = 0, giving Z 2 e 2 Hij(0) — * Sij^j, where 
k 2 — e 2 Z 2 rii/T, the contribution of the i'th ionic species 
to the total Debye screening (wave-number) 2 . We should 
emphasize, however, that (|13fl holds to all orders in the 



2 In the definition 1121 we added an extra factor of T~ 1 , as com- 
pared to the definition in ref. 0] . The development in ref. 0] did 
not assume a uniform neutralizing background, as we did; assum- 
ing instead a neutral system consisting of Boltzmann particles of 
both signs of charge. However the formalism and equations are 
applicable to our case, at the classical level, with the exceptions 
of sum rules based on neutrality of the sea of particles that en- 
ter explicitly. Identical results can be obtained from the usual 
finite temperature QFT approach to the many-body problem; 
we shall explicate the connection in 1201 and the discussion that 
follows. In this latter approach the correlators K(q,u>) and their 
building blocks II are the Fourier coefficients of r-ordered prod- 
ucts (r=imaginary time), and the approximation with the same 
structure as llijl is called the RPA or the ring approximation. 



3 



plasma coupling; a non-perturbative derivation of this 
result is given in rcf. 9] . 3 

In the RPA approximation, we can, of course, employ 
simple graph-summing methods of quantum many-body 
theory, as presented, for example, in ref. 10], to regain 
the classical result l(T3|) . Going beyond the strictly clas- 
sical, we can then add fluctuations in electron density 
to the picture; the degenerate electrons were taken as a 
uniform gas in the derivation of JHJ. At the DH level 
we need only the square of the electron screening wave 
number, 



U\X e \ TV J 



(15) 



where the last approximation is that of complete degen- 
eracy. 

Now we use Ijl3|l for the case of the three species, two 
kinds of ions plus electrons, in the DH approximation, to 
obtain, the ion-ion correlators, 



e 2 q 2 + n\ 



2 2 e 2 q 2 + k\ + k 2 + k\ " [ ' 

We have not written down the expressions for the remain- 
ing three independent elements of the correlator matrix 
K £ye , K lyC i^2,e, since we are still addressing only the ef- 
fects of the electron fluctuations on the baryonic fluctu- 
ations that determine the contribution of the baryonic 
current to the rates. If we were to include the weak cou- 
plings of the neutrinos to electrons, there would be a piece 
of the amplitude from neutrino-electron interactions that 
interferes coherently with the baryonic current, for very 
small values of q, and we would need these correlators 
to calculate the effects. However, we shall argue below 
that values of q for which this occurs are too small for 
the interference to be of interest. 

Using and ljP2l in © we obtain, 



d 7 (q)/dn = {A^GlyEld + cos9)W(q) (17) 



where, 



W(q) = T 



X 2 K 2 (q 2 



e 2 (q 2 +k\ + kI + k 2 ) 



(Ai 



A 2 ) 2 k?k| 



e 2 (q 2 



e 2 (q 2 



(18) 



3 When we go beyond the lowest approximation, the off-diagonal 
parts of the proper polarization matrix do not vanish. 



We note that in the free-particle limit of i|18|) . char- 
acterized by q/n — * oo, and for the case of ions with 
a common value of Z, we then obtain S(q) = 1. But 
q/n remains quite small in any domain of interest in the 
present problem. Note also that we recapture (JSJ as the 
q, n e — > limit of ltTF|l 

We estimate the numerical importance of our of multi- 
component effects for a case in which the density is 10 12 
g c~ 3 , the composition, on the average, is nickel, and 
the temperature is 4 MeV. For algebraic simplicity we 
take the neutron number in both components to be the 
same and also set sinf?^ = 1/2; for the nuclear charges 
we take Z = 28(1 ± S) so that A± = 1 ± 8/2. Then we 
have k 2 w 870 (MeV) 2 and k 2 rj 18 (MeV) 2 . The ratio 
of rate, with our effects included, to the one component 
plasma rate is then given to leading approximation in the 
parameters <5, n e by, 



W(q,5, n e ) 
W(q,0,0) 



1 + [1740 S 2 



18](^) 2 . 



(19) 



Thus when the measure of nuclear diversity, 6, is very 
small, we find a doubling of the one-component plasma 
result for a neutrino momentum transfer of 4 MeV, com- 
ing from the electron density fluctuations. For S = .1 
we obtain a tripling. Larger diversity parameters will 
give much bigger enhancements, as would the choice of a 
smaller value of q. 

The point of this exercise was to evaluate the potential 
impact of our corrections, not to deal with a realistic mix 
of isotopes. In the real problem we have a complex mix 
of nuclei, and there are a wide variety of possibilities for 
the components. In addition, in evaluating the potential 
impact, we need to incorporate the whole scenario of neu- 
trino production through electron capture in the dynamic 
environment just to know which regions of q are most im- 
portant. The part of the spectrum of neutrino momenta 
that dominates the lepton loss rate is dynamically deter- 
mined, and it is clearly centered in a lower energy part of 
the spectrum than the 3 T range that we think about in 
the usual energy transport problems. This comes about 
first through the fact that the opacity is much less for 
the lower energy neutrinos, and secondly because of the 
repopulation of these states through electron capture by 
nuclei. 

Even worse, from the standpoint of using our results, 
is the fact that already in the parameter region which we 
used in our numerical example above, the plasma is mod- 
erately strongly coupled, with a value T = 8.8. In this 
case we expect the D-H results will become inaccurate 
as q is increased beyond a certain point. Determining 
this point requires a computational approach. Looking 
at the "molecular dynamics" results of Luu et al f° r 
the one-component case, plotted in ref. for exactly 
the parameters which we used above, we see, for exam- 
ple, that (for T = 4 Mev) D-H works fairly satisfactorily 
at q = 6 MeV, but is low by a factor of ten at q = 18 
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MeV. Of course, since both the D-H and "molecular dy- 
namics" approaches will give very different results in the 
multi-component case than in the single-component case, 
the guidance provided by the above example could be re- 
garded with suspicion. Thus we believe that an essential 
preliminary to doing real physics in this problem is to 
do Monte Carlo studies of the correlators for a two com- 
ponent plasma, with both the diagonal and off-diagonal 
terms in the density correlators determined with fairly 
high precision, since we now understand that when we 
substitute these results into © there will be near cance- 
lation of the terms, for small q. 

We return to the question of the role of electrons. As 
we saw in the electron density fluctuations do lead 
to some appreciable effects in the case of a single nuclear 
constituent, even in the absence of an electron- neutrino 
coupling. The physics is that the densities of ions and 
electrons fluctuate together. This comes without cost in 
Coulomb energy; this effect is inhibited by the large bulk 
modulus for the degenerate electron sea, but still can be 
significant. 

There are also contributions to opacity from electron- 
neutrino interactions. These have been estimated in 
ref.[l3j. and appear to be numerically small compared 
with the ionic part, in the domain that we used for our 
above estimates. There are two reasons for this: a) the 
ionic terms start with the advantage of the coherence 
factor, of order N (not N 2 , since the electron number is 
larger than the ionic number by approximately a factor 
of N/2); b) electron degeneracy severely limits the final 
electron states available in a neutrino scattering. That 
said, we note that the calculations of ref . |l3( which were 
carried out for an electron gas at zero temperature, prob- 
ably give a significant underestimate of the v-e scattering 
rate. If we are scattering thermal neutrinos from a de- 
generate electron gas at temperature T, there are more 
final electron states available near the Fermi surface, due 
to the diffuseness of this surface, than there would be 
from calculating the volume of the region of non-overlap 
between a T = Fermi sphere and the sphere displaced 
by a momentum, q where q « T '. This becomes even 
more the case for sub-thermal neutrinos. Thus improved 
calculations of the electronic current contributions may 
be in order. 

As we remarked before, we did not include any elec- 
tron scattering contribution in the coherent part of the 
calculation summarized in the result (|18Jl . Strictly speak- 
ing, it belongs there for very small momentum transfers, 
q « T, that is, when the energy transfer is negligible. 
A more analytic form of this remark is that when the 
energy variable is introduced, the multi-component RPA 
equations are still of the form Ijl3|) where Hi t j(q, uS) is 
now energy dependent. However the rates are no longer 
given directly by integrals over Kij(q,u>) but rather by 



integrals involving, 



1 - exp(-w/T) 



lm[K itJ (q,u)] . 



(20) 



When the prefactor in (|20[) can be approximated as 
T/u>, then the energy part of a phase space integration 
gives exactly the integral over the imaginary part of the 
correlator which, through the dispersion relation, pro- 
duces the real part evaluated at zero energy, as in JfjJ. 
But for the values of q that matter in the present case the 
expansion of the prefactor is not justified for the electron 
contributions, because a thermal neutrino colliding with 
a relativistic electron will typically transfer an amount of 
energy of order T. 

In any case, we find that in the regions of interest the 
ionic current and electronic current contributions do not 
interfere very much, at least at the RPA level. We add a 
caveat however; in a strongly coupled plasma we see no 
reason for the rates coming from the two kinds of neutrino 
interaction to separate so neatly. This could provide a 
further complication to a future Monte-Carlo calculation 
aimed at settling the issues raised in the present paper. 

There is a close relation between the above develop- 
ments on neutrino scattering and some important cor- 
rections to Compton scattering in a hydrogen plasma. 
Indeed, the photon-electron interaction that produces 
the Thomson limit is almost identical in form to the 
neutrino-ion interaction in the present paper. The calcu- 
lation again demands the careful consideration of a two 
component plasma, and the mechanics is parallel to that 
presented in the present paper. The effects are actually 
important in the calculation of Compton opacities in the 
solar interior. This subject was discussed in a number of 
references over the years. Boercker 0] carried out the 
calculation that appears to be completely correct, obtain- 
ing significant corrections that are incorporated into the 
solar opacity codes that are in use today. When we go to 
slightly more extreme conditions than those in the solar 
interior, a density of 1 gc -3 and temperature of 10 6 K, for 
example, in a hydrogen plasma, we find, using the analog 
of ©, 



5(<?thcrm) « 5(0) = 



1 



(21) 



where now kj w k e . 4 

Recapitulating some of the conclusions of this paper, 
we found that the q — limit of the structure function is 
generally non-zero, in contrast to the conclusions of the 
large literature that uses "effective" one-component plas- 
mas of some kind. The answer for 5(0), for the case of 



4 In the solar interior we are not in t he r egion in which the limiting 
form q — » can be applied. Ref. Il4l deals with the complete q 
dependence. 



frozen electrons and any number of ionic species, can be 
found simply from energetic arguments using basic sta- 
tistical mechanics, and even when it is relatively small 
it protects against the dramatic suppressions found in 
the current literature. We note that the results are in 
complete contradiction to the results of the procedures 
for ionic mixtures proposed in ref. Q and recently used 
in ref. Q. For application to supernovae, we need the 
extension to finite q in the strong coupling regime. This 
should begin with the numerical investigation of the reli- 
ability of the RPA results i|18|) for the case of a classical 
plasma with two ionic components and for a variety of 
plasma coupling strengths. 
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